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^"v Abstract 

i Ph | In a two-dimensional model of the planetary atmosphere the com- 

pressible convective flow of vorticity represents a strong nonlinearity 
£N . able to drive the fluid toward a quasi-coherent vortical pattern. This is 

«> \ similar to the highly organised motion generated at relaxation in ideal 

I/"") ' Euler fluids. The problem of the atmosphere is however fundamen- 

tally different since now there is an intrinsic length, the Rossby radius. 
f—N , Within the Charney Hasegawa Mima model it has been derived a dif- 

VO [ ferential equation governing the stationary, two-dimensional, highly 

organised vortical flows in the planetary atmosphere. We present re- 

¥2 suits of a numerical study of this differential equation. The most 

• »h j characteristic solution shows a strong similarity with the morphology 

£>^' of a tropical cyclone. Quantitative comparisons are also favorable and 

several relationships can be derived connecting the characteristic phys- 
ical parameters of the tropical cyclone: the radius of the eye-wall, the 
maximum azimuthal velocity and the radial extension of the vortex. 
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1 Introduction 

In the absence of dissipation the two-dimensional models of the planetary 
atmosphere and of the magnetized plasma can be reduced to differential 
equations having the same structure: the Charney equation for the nonlin- 
ear Rossby waves, in the physics of the atmosphere (Charney 1948); and 
the Hasegawa-Mima equation for drift wave turbulence, in plasma physics 
(Hesagawa and Mima 1978). One of the characteristics of these equations 
is the presence of an intrinsic length parameter, which can be seen from the 
fact that there is no space scale invariance, in contrast to the Euler equation 
(see Horton and Hasegawa 1994). This length is the Rossby radius for the 
atmosphere and the ion sonic Larmor radius for the magnetized plasma. 



A long series of observations, experiments and theoretical (analytical and 
numerical) studies has established that the fluids exhibit an intrinsic trend to 
organisation. This is most obvious at relaxation from turbulent states when 
the fluid evolves toward a reduced subset of flow patterns, characterized by 
regular form of the streamfunction as shown for example by Mathaeus et al. 
1991a, Matthaeus et al. 1991b, Kinney et al. 1995, Horton and Hasegawa 
1994 (and references therein). 

It has been shown (Spineanu and Vlad 2005) that the stationary states at- 
tained at late times in the evolution of the Charney-Hasegawa-Mima (CHM) 
fluids are described by the following equation 

Aip H sinh ip (cosh ip — p) = (1) 

zp 

Here ip is the streamfunction and p is a positive constant. This equation has 
been derived within a field-theoretical formulation of the model of interacting 
point-like vortices and reveals that the asymptotic state of organization of the 
physical vorticity in fluids and plasmas is identical to states of self-duality of 
classical field theory of matter interacting with a gauge field. Consistent with 
the physical model from which it is derived, this equation should describe the 
states of fluids/plasmas characterized by the presence of a background of vor- 
ticity (a condensate of vorticity) and a finite intrinsic length (or a velocity of 
propagation of perturbation, like the gravity wave speed or the sonic speed). 
The numerical solutions of this differential equation has provided interesting 
results for the physical systems which it should be able to describe: plasma 
vortices, atmospheric vortices, non-neutral plasma vortex crystals, current 
sheets. Comparisons with experimental results on Navier-Stokes fluid, con- 
sisting of scatterplots (ip,uj) = (streamfunction, vorticity) (de Rooij et al. 
1999), and with the scatterplots obtained in numerical simulations by Seyler 
(1995) are also favorable. A brief discussion of the derivation of Eq.([T]) will 
be presented in Section 2. 

Due to the similar analytical structure (although for largely different mag- 
nitude of parameters) of the plasma vortex and the atmospheric vortex, we 
expect that this equation leads to solutions that may capture the strong 
nonlinear character of the atmospheric vortical flows in those states where 
the stationarity can be assumed a good approximation. Certainly the prob- 
lem of the structure of the atmospheric vortex cannot be reduced in general 
to only fluid nonlinear dynamics, knowing the very important role of the 
heat exchange and moisture transport and condensation processes. These 
are essential elements of tropical cyclogenesis (see Emanuel 1986 and 1989 
and references therein) but it is often accepted that the fluid dynamics is 



an adequate framework to study the vortex structure at the late stage of 
the evolution (Reasor and Montgomery 2001, Kossin and Schubert 2001). 
The fluid dynamic nonlinearity becomes the dominant constraint determin- 
ing the structure of the velocity field when the thermodynamic processes 
have reached the stationary equilibrium. 

The results of a numerical investigation of equation (fTT) , for a range of 
parameters relevant for the physics of the atmosphere, are summarised here. 

The monopolar solutions of the differential equation have the same mor- 
phology as the two-dimensional flow of a tropical cyclone (Willoughby and 
Black 1996, Wang and Wu 2004, Reasor and Montgomery 2001, Kossin and 
Schubert 2001). The solutions are characterized by a very narrow dip of 
the profile of the azimuthal velocity (tangential wind) in the center of the 
vortex. The radius of the "maximum tangential wind" or the radius of the 
eye wall is much smaller than the radius of the vortex. There is a decay 
of the magnitude of the azimuthal velocity toward the periphery, which is 
much slower compared with the fast decay toward the center. We find a very 
low magnitude (almost vanishing) of the vorticity over most of the vortex 
(approx. from the radius of maximum wind to the periphery), while the 
magnitude in a narrow central region is extremely high. Quantitatively, we 
obtain for the diameter of the cyclone's eye a magnitude which compares 
well with the observations. The maximum vorticity is in a realistic range 
and the radial profile of the tangential velocity is similar to what is found 
in observations or with what is obtained in empirical models and numerical 
simulations. These favorable comparisons are valid in the case of mature, 
quasi-stationary tropical cyclones, after the phases of genesis and dynamic 
intensification. 

Besides the structure of the flow, the numerical solutions for a very large 
number of cases provide a large data set from which various correlations 
can be extracted. In this regard we are mainly led to look for correlations 
between parameters that may present interest in practical cases and even- 
tually can be compared with the observations. The main parameters that 
have been collected in the numerical work are: total radial extension of the 
atmospheric vortex, -R max ; the radius where the maximum of the tangential 
velocity is attained, r„ma X ; the magnitude of the maximum tangential wind, 
^max j n addition we have collected the total energy Ef in and vorticity Qfi n 
in the final state flow field. It is encouraging that we find the linear relation- 
ship, Ef in ~ Qfin which has also been revealed in numerical investigations 
of a large ensemble of point-like vortices (Yatsuyanagi et al. 2005). The 
numerical studies carried out until now are organized in the form of several 
scaling relations connecting vortex characteristics to the few parameters of 
the Eq.([T|). They may be used for comparison with observation or with more 



complex theoretical models. 

In addition, the numerical study of Eq. ([T]) reveals the existence of metastable 
states (quasi-solutions) consisting of (a) extremely concentrated vortical flows, 
similar to the cross section of a tornado and (b) collection of vortices with 
symmetric positions in plane (vortex crystals). 

2 The physical model and the field theoreti- 
cal description 

In this section we briefly recall the origin of Eq. ([T]) . It is useful to consider first 
the ideal fluid described by the Euler equation, for which it has been shown 
that it evolves at relaxation toward a very ordered flow pattern, consisting 
of two (positive and negative) vortices. This state persists for long times, 
being limitted by only the effect of some residual dissipation. From numerical 
simulations it has also been inferred the form of the flow function. It has been 
found that the streamfunction ip (x, y) obeys, in these states, the sm/i-Poisson 
equation, an equation which has very special properties. It is an exactly 
integrable equation (by Inverse Scattering Transform) and is connected with 
a wide class of fundamental systems, like the Thirring lattice of spins in plane, 
the affine Toda system, the Gauss-Codazzi equations for embedding a surface 
in three dimensional space, etc. David Montgomery and his collaborators 
have developed a theoretical statistical model which explains the presence 
of this equation in this context (Kraichnan and Montgomery 1980, Fyfe, 
Montgomery and Joyce 1976, Joyce and Montgomery 1973, Montgomery 
and Joyce 1974, Montgomery et al. 1992). 
The Euler equation is 



d 

-^ + (~V ± V x n) • Vj 



w = (2) 



where n is the versor perpendicular on the plane, the velocity and the vor- 
ticity are respectively v = — Vj^ x n and u = r\S7\ip. It is generally 
accepted (but not yet mathematically proved) that the Euler fluid may be 
described by an equivalent model, consisting of a set of discrete point-like 
vortices moving in plane under the effect of mutual interaction. The latter 
is expressed by a potential given by the natural logarithm of the relative 
distance between vortices normalized to the linear extension of the region in 
plane where the motion is bounded (Kraichnan and Montgomery 1980). A 
fundamental observation is that this formulation exhibits a particular struc- 
ture: matter, field, interaction. The matter corresponds to the density of 



point-like vortices in plane; the field is derived from the potantial and can 
be seen as an independent component of the system; the interaction appears 
as the usual combination of the matter current with the potential and leads 
to the equation of motion of the point-like vortices. The sm/i-Poisson equa- 
tion has been derived by formulating the continuum version of the point-like 
vortices model as a field theoretical model of interacting gauge and matter 
fields in the adjoint representation of SU (2) (Spineanu and Vlad 2003). The 
essential point of the latter derivation was the self-duality of the relaxation 
states of the fluid. 

In its simplest form the Charney-Hasegawa-Mima equation is 

JL (i _ vi) i> - [(-V_l^ x n) ■ V.J Vi^ = (3) 

Similar to the Euler equation there is a discrete vortex model for the Charney- 
Hasegawa-Mima equation, where the potential of interaction is now short 
range (the modified Bessel function of zero order), i.e. it decays on the range 
given by the intrinsic length of the CHM equation. This model has been used 
in meteorology by Morikawa (Morikawa 1960) and Stewart (Stewart 1943). 

In a similar approach as for the Euler fluid, it has been developed (Spineanu 
and Vlad 2005) a field theoretical model for the continuous version of the 
point-like vortices with short range interaction, based on the Chern-Simons 
action for the gauge field in interaction with the nonlinear matter field, in 
the adjoint representation of the SU (2) algebra. In this model it is possible 
to derive the energy as a functional that becomes extremum on a subset of 
stationary states and presents particular properties. The general character- 
ization of this family of states is their self-duality, which here means that 
the energy functional becomes minimum when the squared terms in its ex- 
pression are all vanishing, leaving as lower bound for energy a quantity with 
topological meaning, proportional with the total vorticity. 

The result is a set of equations parametrized by the solutions of the 
Laplace equation in two- dimensions. The Eq.flTJ is the simplest of this family. 

3 Numerical studies 

To solve Eq.fJTJ we use the code "GIANT A software package for the numer- 
ical solution of very large systems of highly nonlinear systems" written by 
U. Nowak and L. Weimann (Nowak and Weiman 1990). The code is part 
of the numerical software library CodeLib of the Konrad Zuse Zentrum fur 
Informationstechnik Berlin. The meaning of the abbreviation is: GIANT = 



Global Inexact Affine Invariant Newton Techniques. This code solves non- 
linear problems 

F(V) = (4) 

initial guess of solution, if) = ijj^ ' 

where F (ip) is a nonlinear partial differential operator. The presence of the 
hyperbolic trigonometric functions, with very fast variation with the argu- 
ment, renders the equation difficult to solve and many initial conditions do 
not lead to a solution since they cannot initiate a converging iteration. With 
all the difficulties of finding a right initialization of the integration procedure 
we note however that the solution with the morphology of a tropical cy- 
clone appears insistently from a wide class of initial functions which contains 
vortical flows. 

3.1 Parameters and Initializations 

Eq.(TJE]) has been derived from the self-duality subset of the Euler-Lagrange 
variational equations for the action functional of the field-theoretical model. 
In this action there are only two physical parameters: the coefficient of the 
Chern-Simons action, which we have identified as the sound speed, c s and the 
asymptotic vorticity which is the Coriolis frequency f (or, in plasma physics, 
the ion cyclotron frequency fi c ). The space-like parameter that normalizes 
the Laplace operator in the Eq.flTJ is the ratio 

Pg = C a /f (5) 

i.e. the Rossby radius p g (respectively the sonic Larmor radius in plasma, 
p s = c s /Q C i). All distances implied in the solution are normalized to p g . The 
streamfunction is normalised as 

4' - *£ (6) 

The unit for the streamfunction is p 2 fo and the unit for vorticity /q. Then the 
unit for velocity is p g fo- Here and in the rest of the paper the upperscript phys 
is used to indicate that the quantity is dimensional, measured in appropriate 
physical units. However the units themselves (i.e. p g , /o and combinations) 
are written without this upperscript since there can be no confusion. 

If we know the large radial extension of a vortex (-R^f) the normalized 
parameter -R max is obtained by dividing to p g . The result of integration is 
very sensitive to -R max and this points out the essential role played by p g in 



numerical studies aiming to reproduce observations. When the depth of the 
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atmospheric perturbation is H, the Rossby radius is p g = (gH) ' / f (where 
g is the gravitational acceleration). Actually p g is influenced by several other 
parameters than H, in particular the vorticity and it can have a range from 
tens to thousand kilometers. When combined with the range of vortex ex- 
tensions -R^m, it results that -R max has a range between a fraction of unity 
to few units (in plasma -R max can be hundreads to thousands). For mature 
stationary cyclones the radial extension of the cyclon vortex and the Rossby 
radius are of similar magnitude and this means that we have to study the 
range R max ~ 1. 

The domain of integration is a square with a side of length 2L 

(x,y)e[-L,L]x[-L,L] (7) 

on which we place a rectangular mesh n x n usually with n = 101. 

We impose Dirichlet boundary conditions i.e. the streamfunction is a 
constant, ipj, on the limits of the square of integration. Here ipb is the smaller 
root of the equation cosh?/' — p = and since in all our runs we have p — 1, 
the condition is ipb — 0. This means zero vorticity at the boundary but it 
has defavorable effect on the velocity field, which does not allways vanish at 
the border. 

In general the initial profiles have been of various types: (a) symmetric 
profiles (e.g. Gaussian functions, or various annular shapes) with maximum 
centered on (0,0), or the Petviashvili-Pokhotelov vortex, Eq.fjHj) below; (b) 
functions expressed as product of trigonometric functions; or (c) collections 
of localised vortex-like perturbations placed randomly. In this work we report 
results obtained for initializations with vorticity of the same sign over all the 
space region. It resulted that the monopolar vortex which may be relevant 
for the atmosphere physics can actually be obtained as a solution of the 
convergent iteration from a large class of initial conditions, of various shapes 
belonging to the three classes mentioned above. However, in terms of a 
set of parameters describing any of the classes of initialization, finding the 
convergence (solution) proves to be complicated and finding a solution is 
rather exceptional. It can be said that the regions favorable to convergence 
in the space of initial function are intricate and have sharp limits. 

In order to simulate the formation of a vortex the main series of runs 
reported here have adopted as initial function (from which the iteration of 
GIANT starts) a monopolar, circularly symmetric form characterized by few 
parameters, the Petviashvili-Pokhotelov vortex 

ip (x, y) = ip [sec h (/cr)] 4/3 + ip b . (8) 



Here kr = v/(# 2 + y 2 ) /8 2 , 5 may be seen as a peaking factor of the initial 
shape, ipQ is the amplitude and ipb is the boundary value. The reason to 
choose the form (jHJ) is connected with the set of equation (continuity and 
conservation of momentum) from which the CHM equation is originally de- 
rived. It has been shown by a multiple space-time scale analysis that the 
late stage evolution of the full set of equation is dominated by mesoscopic 
scales (of the order \JL n p s where L n is a typical length of the gradient of 
the equilibrium density) where a different nonlinear mechanism is present. 
Instead of the nonlinearity term in ([3]) there is a scalar or Korteweg-DeVries 
nonlinearity of the type ip-^ in one dimension, and the equation actually is 
replaced by the Flierl-Petviashvili equation (see for example Spineanu et al. 
2004). Or, this equation has a solution expressed as Eq.flHJ). As initial func- 
tion, Eq.© has some advantages: it has physical relevance for the system 
which is behind the stationary states emerging from Eq.([3]); it is a vortical 
structure, as those which may be expected to form spontaneously in real sit- 
uations; and it has few parameters which we take as coordinates in a space 
of initial functions that must be sampled to find solutions. 

The parameters of the equation and of the initial function, which com- 
pletely defines a numerical experiment, are: (a) half the length of the side of 
the square area taken as region of integration, L = L phys /p g ] (b) the constant 
of the equation which in these runs is taken p — 1. This implies that the 
boundary condition for the streamfunction, which also means zero vorticity, 

is ipl = ln f p ± ypM =0; (c) the peaking on the center, described by 

5; (d) the amplitude of the initial function ip . 

The choice of an initial amplitude ipo favorable for reaching a solution 
is made easier if we use the following formula connecting the radius of the 
maximum azimuthal velocity a with the amplitude of the streamfunction at 
this maximum ipo: 

a ~ A/^exp (—0o + 1) (9) 

This formula can be derived by simple manipulations of the equation (JT]), 
together with approximations that are made possible by some numerical 
experience about the orders of magnitude of the normalized quantities. It is 
a rather poor but useful approximation and only works for L < 3, a range 
which is interesting for the atmospheric vortex. 

From the other types of initial profiles we briefly discuss the second one. 
As suggested by the numerical study of the sm/i-Poisson equation the initial 
function has been taken in several runs as a product of trigonometric func- 
tions in both directions, x and y. To have a good initialization, we choose a 
point (x, y) where the initial function is maximum. The equation imposes a 
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condition on only the amplitude ipo of the trigonometric functions of period 
k, 

Aip = tfj [2 (k7i) 2 ] ~ - sinh tp (cosh tp - 1) (10) 

The equation is solved graphically and one of the roots is selected as the am- 
plitude of the initial function. From these initializations we obtain sometimes 
quasi-solutions consisting of multiple vortices. 

3.2 Numerical solution for circularly symmetric vor- 
tices 

When the integration with GIANT identifies a monopolar, circularly sym- 
metric solution (usually with the precision 10 -8 ) the shape is still influenced 
by the square geometry of the domain of integration ((7j) and by the boundary 
condition ipi, = 0. The velocity is usually not zero on the boundaries and this 
means that we have a too approximative value for the total energy (this is 
not the case with the vorticity field and with the energy and vorticity for the 
strongly localized vortices). In these cases we can make a one dimensional 
(radial) integration 

d 2 ip I dip 1 , 

-pr H 7~ + Tr^smhip (cosh-0 — p) = 

dr z r dr 2p z 

with ip = at r = -R max . Since the square and radial problem have different 
boundary conditions the mapping between the solutions obtained in square 
and in radial integrations requires certain care. We have found that the 
general rule 

L = -Rmax/v 2 (11) 



establishes a good correspondence between the square solutions and the radial 
solutions. This simply means that the solution on a square with half side 
L sq = L are very close to the radial solution with extension L rad = i? max 
equal to half the diagonal of the square. This is shown by the Table 1 where 
we compare the important quantities: r„max (the radius where the maximum 
tangential velocity is attained) and t>™ ax (the maximum tangential velocity), 
obtained in square ( "sq" ) and respectively radial ( "rad" ) integrations. These 
quantities are also compared in figures la and lb. 
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Figure 1: Comparison between results on square (squares and continuous 
line) and radial (circles and dashed line) integrations for: (a) v™ 3 * and (b) 
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For radial integrations we use BVPLSQ (Boundary Value Problem Least 
Square Solvers for highly nonlinear equations), written by P. Deulfhard and 
G. Bader. The code is a part of the CodeLib library of Fortran codes sup- 
ported by Konrad Zuse Zentrum fuer Informationtechnik Berlin (ZIB) (see 
the Prolog of BVPLSQ at the CodeLIB web site). Although is more efficient 
than MULCON (also from CodeLIB) or the NAG subroutine D02GAF, we 
have noted that the radial integration has a weaker ability of identification 
of the solution compared with the two-dimensional (square) one, GIANT. 
The solution is obtained at a lower accuracy which we quantify by defining 
a functional error rad 



error 



rad 



d 2 r[\u\-\NL\f 



(12) 



where to is the vorticity and NL is the nonlinear term in Eq.([T]). This is 
shown in a column of the Table. 

A massive series of (automatic) radial integrations has been performed, 
since they are faster than GIANT. Both solutions and quasi-solutions are 
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obtained and it is confirmed that the very precise results of the square inte- 
grations (GIANT) are points on a line of minimum error (fT2|) . They will be 
reported elsewhere. 

4 Results 

4.1 Summary of the numerical results 

The strong nonlinear character of Eq. ([T]) combined with the internal proce- 
dures of GIANT (with no physical significance) imprints a particular struc- 
ture to a space of functions that is explored for exact solutions. The space of 
functions representing initial conditions are divided into disjoint parts, such 
that from one subset one cannot access the final configuration of another 
subset. This means that in order to obtain a particular type of stationary 
(asymptotic) solution one has to initialize in a particular subset, homotopi- 
cally connected to the final state. Most of the initial conditions does not lead 
to convergence and possibly they correspond to turbulent physical states. Ex- 
act solutions are obtained in the form of trivial (i/j = 0) state and monopolar 
or multipolar vortices. 

Besides exact solutions there are sets of functions that are almost solu- 
tions, i.e. velocity field configurations that verify the equation with only low 
accuracy and are normally rejected by the integration procedure at smaller 
tolerance. They are interesting because they appear systematically and ap- 
proximately exhibit the same characteristics for a fixed L. The iteration of 
GIANT gets almost stuck around such a solution, which may suggest that 
they are metastable states of the physical fluid and eventually evolve slowly 
toward an exact solution, a smooth vortex. We call them quasi-solutions and 
we find useful to include them in our discussion. The main reason for accept- 
ing them as interesting and possibly physically relevant structures resides in 
the particularity of our approach: the fundamental object is the action func- 
tional and the configurations described by Eq.([T|) extremize this action. If, by 
indifferently what method (here numerical), it has been possible to identify 
a configuration that seems to be very close to the extremum of the action 
(possibly a local minimum in the function space) then this configuration may 
play a role in the system's evolution. Although a more detailed description 
of the structure of the function space around the extrema of the action is 
still required, the highly concentrated vortices identified numerically seem to 
be close of extremizing the action functional. 

Restricting to the case of monopolar vortices we summarize the results 
by saying that for every L we obtain two types of final vortices: 
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1. smooth, finite amplitude, vortices verifying the equation for any change 
in the accuracy of the integration procedure. There is a unique smooth 
vortex configuration for each L. 

2. quasi-solutions, consisting of strongly localised vortices, with profiles 
for velocity and vorticity that are much narrower compared with the 
smooth vortices. These quasi-solutions are persistently obtained and 
they exhibit a clusterization around a typical shape for a particular L. 

4.2 Solutions with the morphology of tropical cyclones 

We are focusing here on the type of solutions that exhibit strong similarities 
with the morphology of a horizontal cross section of a tropical cyclone. These 
are circularly symmetric solutions with a strong maximum of the tangential 
velocity, strong concentration of vorticity. 

We consider for ilustration the smooth, circularly symmetric vortex ob- 
tained for L = 1.25. Figj2] presents a section of the streamfunction ip (x,y) 
along the diagonal of the domain (x, y). A section of the vorticity u (x, y) is 
presented in FigOl 




Figure 2: Solution ip(x,y) along the diagonal of the square, for the smooth 
vortex at L — 1.25. 



In order to quantify the accuracy of integration we collect in all the do- 
main (x, y) the pairs (?/>, u) and plot them together with the line representing 
the nonlinear term in Eq.flT]), FigJH The scatterplot of {ip,uj) is almost su- 
perposed on this line. The scatterplot of the pairs [a;, —\ sinh ip (cosh ip — 1)] 
(not shown) indicates a close clustering around the diagonal. Other tests 
are possible and they show that the integration is very good on most of the 
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Figure 3: Vorticity calculated from ip(x,y) obtained by integration. 
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Figure 4: Scatterplot (ip,uj), for the smooth vortex at L = 1.25. 
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region and good within the imposed accuracy in the regions where the second 
derivative is very high. 

The tangential component of the velocity is shown in FigJHJ The steep de- 
scent to the center is clearly visible and its radial extension can be compared 
with the extension of the whole domain. 

We have ploted in FigJU] the section along the diagonal of the amplitude 
of the azimuthal component of the velocity. 
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Figure 5: Azimuthal velocity vg(x,y) for the smooth vortex at L — 1.25. 



In general the radial component of the velocity is much smaller than the 
azimuthal component. For this example (L = 1.25) is in a ratio |i> r | / \vg\ ~ 
0.4/4 = 1/10 and integrated over a circle shows no net inflow to the axis. 

The large amount of results for the range of L : < L < 10, allows to 
formulate two remarks. First we note that for larger L the profile of the 
azimuthal velocity shows smaller amplitude and larger radius of the circle of 
maximum velocity. Second, this variation with L is much faster for low L 
(less or comparable to 1). The differences in the quantitative characteristics 
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Figure 6: Azimuthal velocity vg(x,y) along the diagonal of the square of 
integration, for the smooth vortex at L = 1.25. 
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Figure 7: Streamfunction ip(x, y) along the diagonal of the square of integra- 
tion, for the quasi- solution at L = 1.25. 
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Figure 8: Vorticity u(x,y) along the diagonal of the square of integration, 
for the quasi-solution at L — 1.25. 
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Figure 9: Azimuthal velocity vg(x,y) along the diagonal of the square of 
integration, for the quasi-solution at L = 1.25. 

of the vortices for even small variation of L in this range are substantial. 
Below we provide a scaling which shows exponential behavior, Eq.f TTBl . 



4.3 Quasi-solutions : strongly localised vortices 

For comparison we present the profiles of the solution ijj (x, y) FigJTJ vorticity 
u (x, y) FigJH] and azimuthal velocity Vg (x, y) FigJHl for the quasi-solution 
that correspnds to the same L = 1.25. The sections are along the diagonal, 
denoted r. The results for every L's seem to indicate a clusterization of the 
final total energy and final total vorticity. One should note that the total 
vorticity (obtained by integrating over the square) is however smaller than 
that for the corresponding smooth vortex shown in the preceding figures. 

The total energy |p / d 2 r (Vif)) + \ip 2 for a fluid density po, is larger for 
the concentrated vortices compared to the smooth ones. 



4.4 Quasi-solutions : multiple vortices 

There are episodic structures of multiple vortices that are detected as so- 
lutions under a certain precision and which however evlove to symmetric 
monopolar vortex when the system is allowed to run further, under a higher 
precision. 

It is worth to mention that in a numerical experiment we have identified a 
state where two vortices have been formed, placed in symmetrical positions 
along the diagonal of the square domain L = 0.5. The initial function is 
trigonometric with periodicity k = 2 with a coefficient ipo = 3.8. Examining 
this structure with higher precision, after a longer iteration sequence the final 
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solution was again the centered smooth vortex known for L = 0.5. Therefore 
from the point of view of the numerical experience this state of two vortices 
is irrelevant. However, the persistence of this state inside the iterative search 
may indicate that it is close to a solution, possible less structurally stable. 

Four vortices have been obtained in a run starting from trigonometric 
initial function. The initial function is trigonometric with periodicity k = 3. 
The results show the formation of four vortices, as shown by FigJTUl Each 
of them has a structure that is similar to the one presented in FigEl It is 
interesting to note that again the vorticity is almost zero everywhere on the 
domain, except the regions of the four vortices, where it reaches very high 
values. The local tangential velocity presents the same very fast decay to 
the center of the vortex and each vortex is similar in structure with a typical 
cyclone. 
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Figure 10: The scalar streamfunction if)(x,y) for a four-vortices solution. 

5 Scaling laws inferred from numerical re- 
sults 

As mentioned before the nontrivial results of monopolar vortices are sys- 
tematically of two types: a smooth vortex solution and a strongly localised 
quasi-solution. These are always the same for a fixed L. Their characteristics 
strongly depends on L, especially for the lower part of the range, where L is 
few units or less. 
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Figure 11: The contours of the scalar streamfunction if)(x,y) and the vector 



field (v x .v y ) for a four- vortices solution. 



Since the monopolar final states are independent of the initial conditions 
from which we start and of the particular numerical method of solution (GI- 
ANT) the relationships between the characteristics of the final vortices are 
objective and reflect properties of the equation itself. In addition the results 
(smooth vortex and narrow quasi-solution) are unique for a particular L, 
suggesting we can collect all numerical results in a form of nomograms or 
analytic (eventually spline functions) fit. However, we will look instead for 
analytic formulas which, even approximative, are simpler to use. We exam- 
ine (a) the scaling of the maximum tangential velocity with the radius of the 
vortex, Rmax'-, (b) the scaling of the radius of the eye-wall of the atmospheric 
vortex with the length L. Finally we will examine the existence of a linear 
relation between the energy and the vortictym the final states. 



5.1 The relationship between the maximal tangential 
velocity and the extension of the atmospheric vor- 
tex 

We have inferred from numerical data an expression showing variation of 
■u™ ax with the radial extension of the vortex -R max 

V2 



(L) 



a cxp 



/?,, 



(13) 
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for the interval < L = -R max /v^2 < 6. This is shown in FigIT2"l 




Figure 12: The maximum tangential velocity v™ ax as function of L — 
Rmax/V?*- The dashed line represents the fit according to Eq. (I13|) 

For a = 0.97 this scaling works better for L around 1, and is a poor 
approximation over the range less than 0.4. At L > 1.5 this formula over- 
estimates the maximum velocity and will not be used when we dispose, for 
the particular L, of the full numerical set as obtained by GIANT. However 
many observational data fall in the range L ~ 1 where Eq. (113ft is a good fit 
and there it may be useful for a rapid estimation. 



5.2 The relationship between the radius of the eye- wall 
and the extension of the atmospheric vortex 

On the basis of many runs we have tried to infer a possible relationship be- 
tween the radius of the eyewall, (the radius of the circle where the tangential 
velocity is maximum) and the parameter L. 

When the data collected for a larger range of L, < L < 6, is taken 
into account, it appears that there is an approximative linear dependence of 
r v ™a X (L) on L = R mSbK /y/2. 



(14) 



r v rnax (L) = 0.11 ( h L 



As can be seen from FigJT3] this linear fit may present interest for a very 
wide range of L's but it does not work well for low L, < L < 2.5. Or this is 
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Figure 13: The radial position r„max where the tangential velocity v™ ax attains 
its maximum, as function of L. The dashed line represents the linear fit 
according to Eq. iTHl) 



the range that is relevant for the atmospheric vortex. It is necessary to look 
for a different fit for that range. 

First we normalize the eye-wall radius to the radial extension of the vor- 
tex, R m ax ~ s/ili. Then the same numerical information can be organized 
to show the dependence of r v max / R max on the length -R max . The following 
simple function offers a satisfactory fit for the low L range 



Vn.max i 

U B _ _ 

l^m.ax ' 



exp 



Ri, 



(15) 



Although it slightly overestimates the ratio (see Fig fHl) . this formula is prac- 
tical by its simplicity and may be used for estimations based on observational 
data, as will be described later. 



5.3 Energy and vorticity in the final states 

The smooth and respectively the strongly concentrated vortices are persis- 
tently obtained, from a wide variety of initial conditions. This suggests that, 
apart variations due to the accuracy of the iteration process, these states 
really represent a stationary state and respectively a quasi-stationary state 
of the fluid. Their characteristics are completely determined once we have 
fixed the radial extension of the domain -R max or, for square integration, the 
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Figure 14: Ratio of the radius r„m» where maximum of vq is attained to the 
radius of the maximum extension of the vortex, R max represented as function 
of the length L = R mijiX /V2. The dashed line is the fit according to Eq. (IT51) . 

length L, which is half the side of the square on which the integration is 
performed, with if) = on its boundaries. 

In particular the final states for a fixed L are characterised by two quan- 
tities, the total final energy and the total final vorticity and there are two 
pairs (Efi n , Qfi n ): one for the smooth vortex and one for the narrow vortex. 
They are calculated with 

Ephys = ^ c l[R^/V2f (16) 

n x n y 



4 I 

phys ■ 



„,, 2 EE[K^')l 2 + ^>i) 2 ] 



where p$ ys is the density. We denote E the energy E phys normalized to the 
physical coefficient in the first line above. The vorticity is 

Q phys = / (17) 

x ;r;rEE^' J ) 

a 1 = 1 J = l 

When pairs (E, f2)are collected from all numerical experiences for a par- 
ticular L, the precision of numerical determination of the solution is reflected 
in the dispersion of the points around an average one which may be supposed 
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Figure 15: Energy and vorticity in the final states obtained at L — 1.25. 
The dot indicated by the arrow consists of 84 smooth vortices, while the 46 
narrow vortices (open circles) show a dispersion in energy and vorticity. 
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Figure 16: Energy and vorticity in the final states for several low L values. 
The smooth (black squares) and concentrated vortices (small circles) are 
plotted in logarithmic scale. Every plane corresponds to a particluar L value: 
0.25, 0.35, 0.5, 0.8, 1.0, 1.25 and 1.5. 
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Figure 17: (a) Energy and vorticity in the final states corresponding to 
strongly localized vortices. The dashed line shows that there is a linear 
relationship between Ef in and VLf in . The figure (b) expands the region close 
to the origin. 

to be the exact solution. For the smooth vortices the dispersion is insignifi- 
ant, confirming that the numerical identification of this solution is very good, 
for all initializations. For the strongly localized vortex the dispersion is ap- 
parent and is connected with the higher magnitude of the second derivative 
on a very small area. However the results are clearly clustered around an 
average value which one can use to study various scaling relationships. This 
difference is exemplified in Figjl5] showing 130 results for L = 1.25, both 
smooth and narrow vortices, obtained from widely different initial condi- 
tions. The black dot indicated by the arrow actually consists of 84 distinct 
dots representing almost identical results for the smooth vortex at L — 1.25. 
The open circles are 46 quasi-solutions narrow vortices, which have very close 
characteristics {vq (x,y), r^««,etc.). They seem to represent approximations 
of a unique quasi-solution. The dispersion in the final energy and vorticity 
are a consequence of differences in the calculated second derivative of ip on a 
very small area, the reason for which they are actually rejected by GIANT. 
In FigUH] data are plotted for several values of L. 

The first important relation is between the energy E^ n and the vorticity 
Qfin in the final states. The field theoretical model from which the equation 
is derived points out the existence of a lower bound for the energy functional 
for the point-like vortices: the energy is expressed as a sum of squared terms 
plus a supplementary term that has a topological content. The minimum of 
the energy is obtained by taking the squared terms to zero (and this leads 
to the self-dual equations and further to Eq.([T])) and this makes the energy 
equal to the topological term. Integrating over all plane this equality takes 
the form of a proportionality of the total energy and the total vorticity in 
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the fluid motion. Therefore the numerical results collected for all L should 
exhibit a linear relationship between Ef in and Vtf in . Now we have for each L 
two pairs (Efi n , fi/m) an d we can try to verify this linear relations for both. 
The total vorticity can be calculated easily in both cases, since the vorticity 
has a very good spatial limitation around the eye and decays rapidly to zero. 
The total energy Eq. (fToT) includes the integration of the squared velocity 
but the square domain of integration actually does not have everywhere zero 
velocity at the boundary, especially for higher L. Then a certain amount 
of energy cannot be included and Efi n is not reliable. This is the case with 
almost all smooth vortices (except possibly the very small L ~ 0.25. ..0.5, 
where the localization of the smooth vortex is more pronounced). However, 
for the strongly localized vortices this problem does not arise, for any L. 
FigJTT] shows that there is indeed a linear relation between the energy and 
the vorticity. With lower accuracy (which we have mentioned before) the 
radial integrations confirm the proportionality of Efi n and Qfi n - 

5.4 The existence of a threshold in the initial energy 
and vorticity for obtaining a solution 

Numerical simulations of basic fluid equations have shown that there is a 
boundary in the space of the initial configurations which separates two kinds 
of behaviors: on one side there are states from which the fluid evolves to ran- 
dom, turbulent states and on the other side there are initial configurations 
giving in long run organised, highly ordered flow with a vortical pattern. In 
the case of the Navier-Stokes equation this limit has been called the ergodic 
boundary by Deem and Zabusky (1971). In our case there is no time evolu- 
tion and the iterations of GIANT have no particular physical meaning. We 
simply note however that a similar separation occurs in the present case. For 
any fixed L the iteration of GIANT converges to a smooth vortex only if the 
amplitude of the initial function ijj is higher than a particular value, depend- 
ing on the peaking parameter 5. This can be translated into a condition for 
the initial energy and vorticity. Choosing L = 1.25 we have represented in 
FigUHlthe points (E ini , VL in i) corresponding to all types of final results: zero 
(asterisk), smooth vortex (open circles), narrow quasi-solution (black dots). 
For purely orientative purpose a line is drawn, which approximately sep- 
arates initial states leading to trivial solution ip = (at left) from the initial 
states leading to smooth vortices (at right). This is very steep, f2j ni — Vt° ini ~ 
[Eini — Ef ni ) '" with (E® ni , fi° ni ) one of the points on the limit. The other 
line is also orientative, separating the initial conditions leading to smooth (at 
left) and respectively strongly localised (at right) vortices. 
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Figure 18: Set of points in the plane (E ini , Vt irii ) from which the integrations 
for L = 1.25 have been started, with the initial form Eq.([H]). The points are 
distinguished according to the final states: (1) asterisks are initializations 
leading to trivial, zero, final state; (2) open circles are initializations leading 
to smooth vortex solutions; (3) black circles are initializations leading to 
quasi-solutions with strongly localised vorticity. The dotted line is a tentative 
separation of the trivial solutions (at left) from the smooth vortices (at right). 
The dashed line is an orientative separation of smooth from the strongly 
localised vortices. To the right of the righmost black circles there is no 
solution. 
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5.5 Radial profile of the azimuthal velocity: compari- 
son with Holland's model 



In an integration over square domain we obtain (v x , v y ) with very good accu- 
racy but the outer part of the field may be affected by the square geometry. 
For monopolar vortices we repeat the integration in radial geometry, which 
is extended over the length of the diagonal of the square, as explained. This 
profile can be compared with semi-empiric formulas like Holland's or Di- 
Maria. For the Holland's model we use the formula 



1 po V r / 



2 f 2 >. 1/2 



r 2 f 



2 



The parameters are: 6 = 2 (the shape parameter), P c = 990 (hPa) (the 
pressure in the center of the vortex), P e = 1015 (hPa) (the environmental 
pressure), po — 1.15 (kg/m 3 ), f = 5 x 10 -5 (s _1 ). The parameters used in 
this example correspond to 



(r vr ,) phys ~ 10 (km) , R^ = 95 (km) 

and this can be taken as starting point of our calculations based on Eq. (TjQ) 
and the scaling laws derived from it. We calculate the ratio of the two 
distances r£J(l/i2£(^f = 10/95 = 0.1053. We insert this value in the scaling 
described by Eq.(15) and determine the maximum radial extension of the 
vortex (normalised), -R max 



0.1053 



i? max = 1.0931 
or, equivalently L sq = R max /^/2 = 0.7729. The Rossby radius, is 

<£ 95 (km) 
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1 — exp 1 
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-f^rnax 


jD'phys 
-fi-max 



Pg 



Pin 



1.0931 



86.9 (km) 



This is the first physical unit that we need in order to translate our numerical 
results (normalised quantities) into physical quantities. The unit of vorticity 
is /o = 5 x 10~ 5 (s _1 ) and the unit of velocity 

v = pjo = 86.9 x 10 3 x 5 x 10" 5 ~ 4.34 (m/s) 
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The numerical solution of Eq.flT]), obtained with the code BVPLSQ in 
circular symmetry for L rad = -R max = 1.0931 



gives 



9 



This result is confirmed by the integration on the square domain L sq = 
Rm^/V2 = 0.7729, with GIANT, where we obtain v^ x ~ 9.5. As explained 
before, the integration over the square domain is more precise compared 
with the radial one (GIANT versus BVPLSQ). However to compare with 
the profile of Holland's model, we take the result of the radial integration. 
We can now calculate this velocity in physical units 



\phys 



4.34 (m/s) x 9 ~ 39.11 (m/s) 




Figure 19: Azimuthal velocity from the Holland model (continuous line) and 
from Eq.flT]) (circles). 

The Holland's model gives approximately 39.73 (m/s). In addition we 
plot in Fig [19] the radial profiles of the azimuthal velocity from the Holland's 
model (continuous line) and from our Eq.(l) (dotted line), in physical units. 
The similarity of profiles is apparent. We note however that these profiles 
are senzitive to choices that are difficult to measure in observation: for the 
Holland's model the parameters, in particular b, are affected by imprecisions. 
For our model, the main physical input is the radius of the eye wall and 
the maximum extension of the cyclone and these are equally affected by 
imprecisions. 

6 Practical application of the scaling laws 

In this section we discuss how to use the scaling relations we have derived 
when we want to compare with a real observational data of a tropical cyclone. 
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Since for a fixed L the integration provides a unique smooth vortex, all 
data regarding this vortex are available in the form of functions : ip (x,y), 
to(x,y), ve(x,y), and derived quantities : r v max, t>™ ax , Ef in , Qfi n . All are 
normalized and the first task is to identify the physical units that will relate 
these quantities to the physical data. Basically there are two units: p g and 
fo and the latter is assumed known. 

A possible starting point is to estimate (from observations) the ratio 
between the radius of the circle of maximum azimuthal velocity r"x and 
the radius of the minimum disc containing the cyclone, R^ix (here as usual 
the upperscript phys means that the quantities are dimensional). Suppose 
this ratio can be estimated on the basis of a satelite picture. Then we can use 
Eg. (1151) or Fig{T3]to identify the parameter L. Since a real- life observation 
provides also the physical radius of the cyclone, -R^f we can determine 
the Rossby radius from: -R^f = V2L phys = V2p g L where we use L as 
determined. Then the unit of the physical quantities are calculated : p g fo 
for velocity and p 2 fo for the streamfunction. From the runs of the code 
for L we dispose of profiles for all (normalized) variables. Using the units 
we can calculate some characteristics that can be further compared with 
the observations: the maximum azimuthal velocity f^f, the profile of the 
velocity, the vorticity, etc. 

6.1 Example 1 

We use pictures of the profile of the mean tangential wind for the hurri- 
cane Andrew, according to Willoughby and Black (1996). In Fig.3b of this 
reference it is represented the west to east wind profile before the eyewall 
replacement (23 August 1665 UTC). We retrieve the approximate values: 
Vq^ x ~ 68 (m/s), r^max ~ 12 (km) and we assume (with a certain extension 
beyond the limits presented in the figure) -R^f — 120 (km). From this we 
calculate 

phys 

/ „max > 

_^~^± = .1 (19) 

r>phys i on V > 

With this value we turn to Eq. (TT5l) to calculate -R max and then L. 

phys 



max 



D " -nphys 

1 L-max -ftmax 



0.1 (20) 



1 R 

- |1 - exp (^ 

with the result 

R max ~ 1.0217 , L ~ 0.72 (21) 

Taking L = 0.72, we have at this moment at our disposal all the set of results 
that are obtained numerically for the smooth vortex at this L. Coming back 
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to the physical data we now use the spatial extension of the vortex, -R^f — 
120 (km) to calculate the Rossby radius p g , i.e. the space normalization 

= KUl = _J^_ _ U785 , km s 

9 V2L v/2 x 0.72 

Now we can calculate the other physical units: To = fo = 5 x 10~ 5 (s _1 ) (from 
Emanuel 1989, Table 1), for velocity p g fo = 5.9 (m/s), for streamfunction 
P 2 „fo — 0.694x 10 6 (m 2 /s). By numerically solving Eq.([T]) for L = 0.72 we find 
^max _ ;lq g Note that the Eg. (1131) gives a close value: 10.67. The physical 
value results t>™ ax = 64.31 (m/s). This is comparable to 68 (m/s) the value 
shown by Willoughby and Black. The maximum vorticity is |u; max | = 491.2 x 
f = 0.0246 (s- 1 ). 

6.2 Example 2 

For this example we adopt the following input data: the ratio of the full 

spatial extension of the hurricane to the radius of maximum tangential wind 

is -R^f/r^max ~ 9; and the physical extension of the hurricane is -R^f ~ 

300 (km). The data can be compared with the picture taken by NASA at 28 

August 2005, when the hurricane Katrina was above the Mexic gulf, but the 

identification of physical data is certainly approximative. Eq.t fTBl) is used to 

find 

-, / r P h v s \ 

I / ' „max \ 

^ti^'H 1 - 4 ^ - 8313 

It results the Rossby radius p g = H&g/ (V2L) ~ 212 (km). The unit of 
vorticity is the Coriolis parameter To = f = 5 x 10~ 5 (s _1 ) and we have 
the unit of velocity v = p g To = 10.6 (m/s) . Looking again to the results 
from the numerical integration for L = 0.83 we find the magnitude of the 
normalized tangential velocity v™ ax = 8.37 (note that Eq. (ll3l) gives a similar 
value, 8.26), which means that in physical units we have 

<T X ~ 88.6 (m/s) 

This gives a very high value for the maximum tangential wind, but the 
range is still realistic. The maximum vorticity is |^ max | — 290.38 x fo — 
0.0145 (s- 1 ). 

6.3 Example 3 

We take for the third example the approximative value R^^ /r^Jii ~ 8 and 
a radius of maximum extension of the hurricane of -R^f ~ 350 (km). (This 
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is inspired by the picture taken by NASA on the hurricane Rita, September 
21, 2005). From this data it is obtained 

, / r phys \ 

I / ' „max \ 

L = (-2) In 1- 4^— = 0.9803 

V2 \ I&J 

Then we can calculate the Rossby radius p g = R^ax/ (V^L) = 252.47 (km). 
Taking the Coriolis parameter /o = 5x 10~ 5 (s^ 1 ) we have the physical unit 
of velocity v = p g ZJ = 12.6 (m/s) . The maximum tangential velocity for 
L = 0.98 can be calculated using the scaling formula Eq. (fT3l . with the result 
^max _ g 24. We will use however the exact result, provided by numerical 
integration with GIANT of Eq.(pQ), i>™ ax = 6.15 and this leads to the physical 
velocity 

V 0mL = 77 - 5 ( m / S ) 

The maximum vorticity is \u> max \ = 158.34 x f = 0.0079 (s _1 ). 

These examples confirm the possibility of using the scaling relationships 
Eq. tTTBl and Eq. (TT5l) to estimate the physical charactersistics of the tropical 
cyclone. We note that the results of the estimations can be significantly 
affected by the approximations on the observational data, especially that of 
the ratio r„max / i? max . This is because the magnitude of L which is obtained 
from Eq. ffTol) using a reasonable input value for this ratio belongs to a range 
where the variations with L of all the characteristics of the solution (if), u, 
vg, r„max) are substantial. It is sufficient to look at the dependence of w™ ax 
on L, 9 FigfIU 

7 Discussion 

At the origin of our approach it is the Charney-Hasegawa-Mima model, a two- 
dimensional, nondissipative and purely fluid-dynamical (no thermal process) 
model. Although is a simplified model it exhibits (via the field-theoretical 
formulation) a compact analytic and algebraic structure, self-duality, leading 
to Eq.(0Q). We have several arguments in favor of the conclusion that Eq.(CE]) 
may represent the fluid nonlinear-dynamic part of the atmospheric vortex. 
First, the profiles obtained by solving Eq.(pQ) are similar to results already 
known (from observations or numerical simulation) for the same quantities. 

1. The profile of the the tangential velocity, our Fig El is similar to typical 
tropical cyclone velocity profiles, as represented in Fig. 2 of Wang and 
Wu 2004. This is also confirmed by close similarity with the Fig. la from 
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Reasor and Montgomery 2001 and with the profile given by the Holland 
model or with the experimental observation (for Andrew hurricane, 
Willoughby and Black 1996). 

2. The fast decay to zero of the vorticity u shown in our FigfSjis similar to 
what is shown in Fig. la of Kossin and Schubert 2001; ideally (without 
dissipation and friction) for mature hurricanes the maximum of oj is on 
the center. 

3. We note that in a series of reported numerical simulations, the ten- 
dency of the fields is to evolve toward profiles that are very close to 
those shown in our figures [2] and [5j For example, the Fig. 7a and 7b 
of Kossin and Schubert 2001 show the evolution of the vorticity and 
mean of the tangential velocity from initial profiles which correspond to 
a narrow ring of vorticity to profiles that show clear ressemblance with 
our figures [2] and El The same striking evolution to profiles similar to 
ours appears in Figs. 7 a and b of the same Reference. We have investi- 
gated whether a radially annular profile of vorticity can be a solution of 
our equation ([T}). The result is negative, which may explain why such 
an initial profile evolves to either a set of vortices (vortex-crystal) or 
to a centrally peaked structure as in Figj2j 

Second, we obtain a good consistency between our quantitative results 
for an atmospheric vortex (using approximative input information) and the 
values measured or obtained in numerical simulations. 

Finally, the quasi-solutions which appear to be an interesting feature of 
this equation, are compatible with a series of previously known results: the 
four vortices represented in our figure [TO] are similar to the Figure 4a from 
the work of Kossin and Schubert 2001. And the evolution to a monopolar 
structure we obtain is similar to the same process reported in this reference. 

The numerical results have made possible to formulate several scaling 
relationships connecting parameters of the atmospheric vortex. These are 
simple formulas, intended for practical use, inevitably approximative. This is 
because our purpose was to infer an analytical expression for a curve which is 
determined numerically (or: is known in the form of a simple table of values). 
Other expressions are possible and are worth to look for. A true scaling law 
will become possible only as a result of an analytical investigation of the 
properties of the equation. On the other hand, this equation does not have a 
Backlund transform and is probable not integrable by the Inverse Scattering 
Transform. 
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An important field of future investigation (analytical and numerical) is 
represented by the quasi- solutions, elements of a space of functions that are 
very close of verifying the equation. Are they close to the extremum of the 
action functional of the full field-theoretical model, for example, are they 
metastable states? From the physical point of view it may result that these 
states can be rendered more stable by processes that are connected with 
what is missing from the Charney-Hasegawa-Mima model: third dimension, 
viscosity, thermal processes. It is worth to examine the role of the strongly 
localized quasi-solutions in models of tornadoes. 

An interesting suggestion results from the massive series of radial integra- 
tions: it appears that in the space of functions there are strings of vortical 
quasi-solutions emanating from the exact solution and showing increasing de- 
gree of departure from exactness (increase of the error functional). If some 
physical factor will provide stability to these quasi-solutions, then the sys- 
tem may slide along this string, with the consequence that there is no pure 
stationarity but a continuous evolution to stronger and stronger localisation 
of the vortex. This requires further study. 

Developing from the present one, a future self-consistent model will have 
to include variation of the Rossby radius with the dynamical properties of 
the vortex. Since this implies to consider that the coefficient of the Chern- 
Simons part in the Lagrangian is a nonlinear function of the scalar field, it 
is difficult to say if the self-duality will be maintained. 

The investigation of this equation, and, most important, of the field- 
theoretical model from which it is derived, are worth to be continued. 
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